Spin dynamics in the diluted ferromagnetic Kondo lattice model 
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The interplay of disorder and competing interactions is investigated in the carrier-induced ferro- 
magnetic state of the Kondo lattice model within a numerical finite-size study in which disorder is 
treated exactly. Competition between impurity spin couplings, stability of the ferromagnetic state, 
and magnetic transition temperature are quantitatively investigated in terms of magnon properties 
for different models including dilution, disorder, and weakly-coupled spins. A strong optimization 
is obtained for T c at hole doping p « x, highlighting the importance of compensation in di- 
luted magnetic semiconductors. The estimated T c is in good agreement with experimental results 
for Gai_ x Mn x As for corresponding impurity concentration, hole bandwidth, and compensation. 
Finite-temperature spin dynamics is quantitatively studied within a locally self-consistent magnon 
renormalization scheme, which yields a substantial enhancement in T c due to spin clustering, and 
highlights the nearly-paramagnetic spin dynamics of weakly-coupled spins. The large enhancement 
in density of low-energy magnetic excitations due to disorder and competing interactions results 
in a strong thermal decay of magnetization, which fits well with the Bloch form Mo(l — BT 3 ^ 2 ) 
at low temperature, with B of same order of magnitude as obtained in recent squid magnetization 
measurements on Gai_ x Mn x As samples. 



I. INTRODUCTION 

The discovery of ferromagnetism in diluted magnetic 
semiconductors (DMS) such as Gai_ x Mn x Asji*2i with 
transition temperature T c ~ 110 K for Mn concentra- 
tion x ~ and ~ 150 K in films for x in the range 
6.7 — 8.5%^ has generated tremendous interest not only 
in view of potential technological applications, but also 
due to the novel ferromagnetism exhibited by these sys- 
tems in which magnetic interactions between localized 
spins are mediated by doped carriers. The long-range 
oscillatory nature of the carrier-mediated spin couplings 
results in a variety of interesting behaviour such as sig- 
nificant sensitivity of spin stiffness and transition tem- 
perature T c on carrier concentration, competing antifer- 
romagnetic interaction and noncollinear ordering, spin- 
glass behaviour, spin clustering and disorder-induced lo- 
calization etc.iVAVVMVVVVe,! WV" ag recently 
reviewed iSLSS 

DMS such as Gai_ x Mn x As are mixed spin-fermion sys- 
tems in which the 5 = 5/2 Mn ++ impurities replace 
Ga +++ , thereby contributing a hole to the semiconduc- 
tor valence band. However, large compensation due to 
As antisite defects reduces the hole density p to nearly 
10% of Mn concentration x, which plays a key role in the 
stabilization of long-range ferromagnetic order, and also 
provides a complimentary limit to Kondo systems. The 
interplay between itinerant carriers in a partially filled 
band and the localized moments is conventionally stud- 
ied within a diluted ferromagnetic Kondo lattice model 
(FKLM), wherein — JSi.tr j represents the exchange in- 
teraction between the localized magnetic impurity spin 
Si and the itinerant electron spin <T/. 

Recently, finite-temperature spin dynamics due to 
thermal spin-wave excitations has been studied in 
Gai-a-Mn^As samples with different Mn content (thick- 



ness about 50nm, Mn content ranging from 2% to 
6%) using SQUID (superconducting quantum interfer- 
ence device) magnetization measurements^ 3 - The tem- 
perature dependence of (low-field) spontaneous mag- 
netization shows nearly linear fall off, similar to ear- 
lier results exhibiting even a distinct concave behaviour 
for unannealed samples* 3 ^ possibly resulting from spin 
re-orientation transitions due to temperature-dependent 
magnetic anisotropics^ However, the spontaneous mag- 
netization obtained using linear extrapolation from a 0.3- 
0.4 T magnetic field to overcome the anisotropy fields, 
discussed earlier for epitaxial ultra-thin Fe and FeCo 
filmsj2£ is found to be well described by the Bloch form 
M(T) = M (l — _BT 3 / 2 ), with a spin- wave parameter 
B ~ 1 — 3 x 10~ 3 K~ 2 / 3 which is about two orders of mag- 
nitude higher than for Fe and FeCo films. This large dif- 
ference cannot be attributed only to a reduced exchange 
interaction. Post-growth annealing has been shown to 
significantly increase T c , possibly due to enhancement of 
carrier concentration resulting from decrease of Mn inter- 
stitial concentration^ For a 50nm sample with 6% Mn, a 
decrease in B from 2.7 x 10" 3 K" 2 / 3 to 1.4 x 10" 3 K" 2 / 3 
has also been obtained upon annealing^ with a corre- 
sponding increase in the spin-wave stiffness constant D 
from 53 meVA 2 to 71 meVA 2 , which is of same order 
of magnitude as obtained from magnetic Kerr measure- 
ments using pump-probe setup of standing spin waves in 
ferromagnetic Gai_ x Mn x As thin filmsi^ 

In view of these recent findings of strong thermal de- 
cay of magnetization in DMS systems, it is of interest 
to investigate the interplay of disorder and competing 
interactions on magnon excitations which typically con- 
stitute the lowest-energy excitations and therefore essen- 
tially determine the low-temperature behaviour in mag- 
netic materials. In this paper we consider the ferromag- 
netic Kondo lattice model (FKLM) with several exten- 
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sions; we treat disorder exactly by considering explicit 
disorder realizations on finite-size systems (random im- 
purity locations, random on-site potential etc.) and av- 
erage over sufficiently large number of configurations to 
obtain statistically reliable results. We also quantita- 
tively investigate finite-temperature spin dynamics due 
to thermal excitation of magnons within a locally self- 
consistent magnon rcnormalization scheme, equivalent 
to a site-dependent Tyablikov decoupling (local RPA), 
and present the first site-dependent calculations for lo- 
cal impurity magnetization by explicitly incorporating 
the spatial feature of magnon states. As we shall see, 
the disorder-induced formation of low- and high-energy 
localized magnon modes, corresponding to weakly- and 
strongly-coupled spins respectively, results in a variety 
of interesting spin-dynamics behaviour, such as con- 
cave magnetization behaviour due to dominant nearly- 
paramagnetic contribution of weakly-coupled spins and 
an enhancement in T c due to strong local correlations in 
impurity-spin clusters. A simplified calculation involv- 
ing two discrete spin-excitation energy scales and a finite 
fraction of weakly-coupled spins was discussed earlier ^£ 

Magnon excitations provide a composite measure of 
the carrier-induced spin couplings in the collinear ferro- 
magnetic state, with negative-energy modes signalling in- 
stability due to competing antiferromagnetic (AF) spin 
interactions. Magnon properties have been studied as 
function of electron density n in the conduction band 
and the spin-fcrmion coupling J within the concen- 
trated FKLM (having a magnetic impurity at every lat- 
tice site) in the context of heavy fermion materials (21 
ferromagnetic metals Gd, Tb, Dy, doped EuX2& and 
manganitesiS^Si^iiS In the context of DMS, magnon 
properties have been studied earlier in the random phase 
approximation (RPA) for the impurity-band modeUS 
and for the diluted Hubbard modeli&Sii where disor- 
der was treated exactly within finite-size numerical stud- 
ies, and for the diluted FKLM within the virtual crys- 
tal approximation (VCA) where a uniform impurity- 
induced Zeeman splitting of the carrier spin bands 
is assumed^ within the coherent potential approxi- 
mation (CPA)jiSiS and also for ordered impurity ar- 
rangements to make quantitative comparisons with dif- 
ferent approximations^ Magnon spectrum and tran- 
sition temperature have also been obtained recently 
for Gai_ x Mn x As and Gai_ x Mn x N in terms of effec- 
tive Heisenberg models with realistic exchange cou- 
plings, obtained recently from first-principle calculations 
as welli 35 i 36 

Within mean-field (MF) approaches for the magne- 
tization behaviour in DMSr2*iSii& it is the large sep- 
aration of spin and fermion energy scales which plays 
the key role. For hole concentration p in the concen- 
trated FKLM, the very low effective ordering field for 
spins ( Jp) compared to that for fermions (JS) results in 
an intermediate-temperature regime Jp <C k^T <C JS 
wherein the nearly-paramagnetic contribution of spins 
accounts for the unusual concave shape ,i2ii£ii2i Every im- 



purity spin contributes equally in this picture and disor- 
der has no direct role, as in other VCA-based magneti- 
zation studies within spin-wave theoryi2ii£ However, the 
magnon energy scale is typically an order of magnitude 
smaller than the MF low-energy scale Jp, even in the dis- 
ordered case where fermion spin polarization is substan- 
tially reduced in impurity-poor regions, implying that the 
relevant energy/temperature scales for spin dynamics are 
actually determined by these collective excitations. Un- 
like the purely local mean field Jp, the collective excita- 
tions involve non-local spin couplings Jy , and hence not 
only incorporate competing interactions and frustration 
effects, but also differentiate between strongly-coupled 
cluster spins and weakly-coupled isolated spins. As only 
the finite fraction of weakly-coupled spins involved in the 
low-energy localized magnon states yield the nearly para- 
magnetic contribution, this picture is quite different from 
the MFT picture. 

The organization of this paper is as follows. The RPA- 
level theory for magnon excitations in real space is de- 
rived in section II for a general fermion Hamiltonian, in- 
cluding dilution, disorder, and multiple bands, and the 
Goldstone-mode behaviour expected from spin-rotation 
symmetry is explicitly verified. Results for several dif- 
ferent ferromagnetic Kondo lattice models — involving 
dilution, disorder, and weakly-coupled spins — are then 
discussed in Sections III, V, and VI, respectively, with 
finite-temperature spin dynamics introduced in section 
IV. Conclusions are presented in section VII. 



II. MAGNON EXCITATIONS 

Magnons represent transverse spin fluctuations about 
the spontaneously broken-symmetry state and consti- 
tute gapless, low-energy excitations for magnetic sys- 
tems possessing continuous spin-rotational symmetry. At 
low temperature, magnons therefore play an important 
role in diverse macroscopic properties such as existence 
of long-range order, magnitude and temperature depen- 
dence of the order parameter, magnetic transition tem- 
perature, spin correlations etc. In the following we con- 
sider finite temperature T, and obtain magnon excita- 
tions at the RPA level where magnons interactions are 
neglected. Finite impurity concentration, disorder, and 
interaction strength are treated on an equal footing in 
this approach. 

We consider the Kondo lattice model 

Jf^o-^V, (1) 
i 

where Hq represents the free-fcrmion part consisting gen- 
erally of hopping and on-site energy terms, and the sec- 
ond term represents the exchange coupling between im- 
purity spins Sj and fermion spins <Tj/2 at impurity sites 
/. The analysis for magnon excitations discussed below 
is independent of details of Hq, for which several differ- 
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ent cases including hopping and potential disorder are 
considered in later sections. 

Applying the approximate Holstcin-Primakoff trans- 
formation from the spin- lowering (Sf) and spin-raising 
(5/) operators to boson (magnon) creation and annihi- 
lation operators b\ and bj, 



Sy = b\y/2Sl 
Sf = Sj-b^b! 

the Kondo lattice Hamiltonian reduces to 



H = H 



blCTj 



(2) 



(3) 



where erf = erf ± ier V j and the "spin quantum num- 
bers" Si = (Sj)mf refer to finite-temperature magneti- 
zations obtained self consistently in the mean-field state. 
The above approximate transformation neglects quartic 
magnon interaction terms of order 1/5. 

Starting with a MF approximation, (b\) = (bj) = 

(bjbi) = 0, the Hamiltonian (3) decouples into a fermion 
part with an impurity-field term 



^fermion — ^0 



J 



and a local boson part 



•x/0 
''boson 



I 



b% 



(4) 



(5) 



representing the energy cost of a local spin deviation. 
At the MF level, determination of impurity and fermion 
magnetizations (Sf) and (<r|) involves a self-consistent 
solution of the coupled spin-fcrmion problem in terms of 
Brillouin and Fermi functions. In the zero-temperature 
limit, as Si = (S^ )mf = 5, the impurity magnetic field 
seen by fermions has same magnitude JS/2 on all sites, 
however, the magnetic field J(ct|)mf/2 seen by impurity 
spins remains non-uniform due to positional disorder. 

Proceeding next to transverse spin fluctuations about 
the MF state, we obtain the time-ordered magnon propa- 
gator for the impurity spins in terms of the corresponding 
boson propagator 



g Ij -(t-t) = i(* G \T[sj(t)Sj(t')] 

= v/257 (»<* G | TMttfjtf)] 




(6) 



at the RPA level by summing over all bubble diagrams 





FIG. 1: Distribution of fermion spin densities on impurity 
sites, showing significantly reduced densities in impurity-poor 
regions relative to their average value « p/x. Here J = 4, hole 
doping p » 4% and impurity concentration x » 20%. 
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with Fermi functions /jf and f m i involves integrating 
out the fermions (eigenvalues {Ei a } and wave functions 
{ipia}) in the broken-symmetry state. It is the particle- 
hole bubble [x°(a>)]/j which mediates the carrier-induced 
impurity spin couplings in the ferromagnetic state, and 
the oscillatory, long-range nature of the spin couplings is 
effectively controlled by the fermion band filling and the 
impurity field strength. 

In terms of the site-diagonal zeroth-order magnon 
propagator 



= 



[25/ 



boson 
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the full magnon propagator can then be expressed as 



i 



£ [x°M][g°M] 
i 



w-[tt(w)L 
in terms of a boson "Hamiltonian" 

[H(u))]u = SiSu - Su(u) 
involving the boson on-site energy 
J , 



£ 



(9) 
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and the boson hopping terms 
'J 2 



■[xV)]i 



(12) 



associated with carrier-induced spin couplings Jjj ~ 
(J 2 /4)[x°(u})]ij ■ While dynamical effects are in principle 
included in the magnon Hamiltonian [7i(w)], we find that 
the to dependence is sufficiently weak to be neglected, 
so that the eigenvalues and eigenvectors of \H\ directly 
yield the (bare) magnon energies {w ; } and wave func- 
tions {0°}. 

Equation (9) has exactly same structure as obtained 
from the Tyablikov decoupling for an effective Heisen- 
berg model with spin couplings Ju ~ ( J 2 /4)[x°(w)]/j, 
and readily yields a locally self-consistent renormalized 
magnon theory, as discussed in section IV. 

In order to obtain a zero-energy Goldstone mode (bo- 
son amplitude ^/2Si) consistent with spin-rotation sym- 
metry, the energy cost of creating a local spin deviation 
must be exactly offset by the delocalization-induced en- 
ergy gain, which requires that 



which is indeed exactly satisfied. This is easily verified 
in the concentrated limit where translational symmetry 
results in plane-wave fermion states with band energies 
Eka = £ k — aJ(S z )/2, and the particle-hole propagator 
(7) simplifies to 

X °(q,u = 0) = ^[xV = 0)]/j 



= E 
k 

= E 



/ kT (i-/ki) , (i-/ kT )/ki 



J(S Z 



-J(S Z 



(/k T ~ /ki) _ (Q 
J(S Z ) J(S Z ) 



(14) 



which ensures that the required condition (13) is satis- 
fied. The required condition (13) for a spin-rotationally- 
invariant system can be easily derived from a pertur- 
bation analysis for the transverse fermion-spin density 
due to small transverse impurity fields (corresponding to 
small twist of the spin coordinate system) and then com- 
paring with that expected on symmetry grounds. 

The transverse spin-fluctuation propagator for the 
fermion spin, defined in analogy with (6) in terms of 
fermion spin operators a~ and tr + , involves the same 
bubble sum at the RPA level and is given by 



[x»] 



1 + £[0°(«)][X»] 



(15) 



While the fermion spin-fluctuation propagator also in- 
volves single-particle (Stoner) excitations, described by 
the poles of [x°(w)], the collective (magnon) excita- 
tions are, as expected, same as for the impurity spin- 
fluctuation propagator, as it is the fermions which medi- 
ate the impurity spin couplings. We note that the above 



fermion propagator involves an effective local (but dy- 
namical) fermion interaction 



[u{u>)] = ^[g\u)\ 



(16) 



corresponding to the local magnon exchange. Fur- 
thermore, an effective isotropic fermion interaction 
U(uS)(Ti.<ji/A, expected from spin-rotational invariancc 
in the fermion sector, provides not only the required 
transverse terms (15,16) but also the correct HF descrip- 
tion (4) from the longitudinal term crjaj. These obser- 
vations allow for a spin-rotationally-symmetric analysis 
of quantum corrections to the magnon propagator due 
to self-energy and vertex corrections along same lines as 
studied recently within the Hubbard models 

For a multiband case, with an exchange coupling 
— i ^2 Ia J a Sj.crj a involving sum over the fermion band 
index a, a straightforward generalization of the above 
analysis yields the same structure (9) for the magnon 
propagator, with magnon on-site and hopping energy 
terms given by 



a l)a 



(17) 



(18) 



a/3 



If the free-fermion Hamiltonian Ho does not involve any 
band mixing, as considered in recent studies^ then 
[X°(w)]"j is diagonal in band indices, and [x°( a; )]?j f° r 
band a is simply given by Eq. (7) in terms of states for 
the corresponding fermion band. A recent equation-of- 
motion analysis for the multiband case yields a similar 
structure for the mRKKY spin couplings^ 

The dimension of the magnon Hamiltonian Ti is N m , 
the number of magnetic impurities per configuration. 
From the N m magnon energies {to®} and wave functions 
{$} we evaluate the magnon density of states and par- 
ticiptation ratio (PR) 



1 1 



N(lu) = 

7T iV m 

PR = 1/J2 



E 



V 



(ijj — LOj)' 2 + 7] 



(v->o) 
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h ) 



(19) 



which together provide a complete picture of both spec- 
tral and spatial features of magnon states. The participa- 
tion ratio provides a measure of the number of impurity 
spins over which the magnon state is extended. Gener- 
ally, if the normalized wave function 4>i corresponds to a 
state with essentially non-zero amplitude (f>\ ~ l/\/n on 
n sites, then PR~ n. The participation ratio for magnon 
states can thus range between A^ m for a fully extended 
magnon state and 1 for a site-localized magnon state. 
The PR thus readily allows localized magnon states 
(PR~ 1) to be distinguished from extended magnon 
states (PR~ N m ). 
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In order to quantitatively study the interplay of disor- 
der and competing interactions, and the role of disorder 
on spin dynamics we consider below different disorder 
models for the fermion Hamiltonian Hq on finite-size lat- 
tices with explicit disorder realizations of impurity po- 
sition, on-site energy etc. so that disorder effects are 
treated exactly here. We consider throughout a simple 
cubic lattice for simplicity, with periodic boundary con- 
ditions. The number of host lattice sites N — L 3 de- 
termines the dimension of the fermion Hamiltonian H 
to be diagonalized, and we have considered system sizes 
L = 8, 10, 12. In all cases we consider the saturated ferro- 
magnetic state with number of spin-"f fermions N<\ = N 
and doping in the spin-| band which is pushed up by the 
impurity Zeeman field. The impurity and carrier (hole) 
concentrations referred below correspond to x = N m /N 
and p = (N — N±)/N, respectively. 

III. DILUTED KONDO LATTICE MODEL 

Providing a minimal description of the exchange cou- 
pling in DMS systems such as Gai_ x Mn x As between 
Mn impurity and carrier spins, the diluted Kondo lat- 
tice model 

H = t a\,a a i+s,o- + £d ^ a \,<r a i,° ~ 2 S/ - cr/ ( 20 ) 

J , (5.CT 7,fT I 

represents N m magnetic impurities placed randomly on 
a fraction (/) of the N host sites (i), with impurity con- 
centration x = N m /N. We consider a positive nearest- 
neighbour hopping t so that the host k = state lies 
at the top of the valence band; doped carriers (holes) 
go in long wavelength states, so that small-fc particle- 
hole processes near the Fermi energy are dominant in the 
carrier-induced ferromagnetic spin couplings, and there- 
fore other details of the energy band are relatively unim- 
portant. In the following we set t = 1 as the unit of 
energy scale. Also, as the temperature regime of inter- 
est, set by the magnon energy scale, is very low com- 
pared to the MF energy scale, in the following we have 
only considered the T = case which provides a good 
approximation of the low-temperature MF state. 

Fermions see an effective potential disorder due to di- 
lution; both disorder and dilution cannot be changed in- 
dependently. We have therefore included an impurity 
on-site energy (chemical potential) e^, which provides an 
effective control of disorder independently of dilution. For 
negative e^, the reduced impurity-site mean-field energy 
(e<j + JS/2) for spin- J, fermions reduces the disparity be- 
tween host and impurity sites, resulting in lower effective 
disorder (for = —JS/2, there is no disorder at MF 
level!). The effects are quite dramatic on the magnon 
spectrum, showing significant decrease in the density of 
low-energy modes and hence enhanced stability of the 
carrier-mediated ferromagnetic state. 

The fermion spin polarization (af) typically shows sig- 
nificant site variation as fermions tend to accumulate in 
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FIG. 2: Magnon density of states for different impurity dilu- 
tions for a N = 8 3 system with J — 4. The upper and lower 
panels correspond to hole doping p m 4% and « 11%, respec- 
tively. Inset shows (for x ~ 20%) comparison of magnon DOS 
for different system sizes N = L 3 . 



impurity-rich regions, leaving significantly reduced polar- 
ization in impurity-poor regions. Figure 1 shows the dis- 
tribution of fermion spin polarization on impurity sites, 
obtained by diagonalizing the fermion Hamiltonian (4) 
on a N = 8 3 system for 50 configurations. Besides the 
broad peak at the expected value ofp/x corresponding to 
average hole density per impurity site, there is an addi- 
tional peak at significantly reduced fermion spin polariza- 
tion corresponding to impurity-poor regions; this intro- 
duces a new low-energy MF scale J(af)S/2 such that at 
comparable or higher temperatures these weakly bound 
impurity spins become nearly paramagnetic. However, 
as we shall see, the magnon energy scale is an order of 
magnitude smaller than even this low-energy MF scale, 
indicating that the dominant spin dynamics is due to 
thermal excitation of the collective magnetic excitations 
rather than that described by the Brillouin function cor- 
responding to single-spin excitation energies within the 
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FIG. 3: Participation ratio for different impurity dilutions 
for a N = 8 3 system with J = 4 and hole concentrations 
p ~ 4% and ~ 11% in upper (four) and lower (four) panels, 
respectively. 



FIG. 4: Magnon density of states for different J and td, show- 
ing effective control of impurity disorder independently of di- 
lution, for a N — 8 3 system with p ~ 4%, x ~ 40%, and 
€d = and J = 4.0 in upper and lower panels, respectively. 



MF theory. As the temperature regime of interest is 
therefore much lower than this MF energy scale, a self- 
consistent finite-temperature mean-field analysis is not 
required and therefore the T = description provides a 
good approximation of the low-temperature MF state. 

Figure 2 shows the configuration-averaged magnon 
density of states (DOS) for two different hole doping con- 
centrations and for varying degree of dilution. Initially 
dilution is seen to essentially broaden the magnon spec- 
trum, but higher dilution results in significant softening 
of the collective excitations, with the magnon DOS peak 
progressively shifting to lower energy. Furthermore, this 
magnon softening is distinctly more pronounced at lower 
hole doping, reflecting more effective competition due to 
longer range of spin couplings induced at lower hole dop- 
ing. In addition, new structure appears at higher energy, 
particularly at higher dilution, which is due to localized 
magnon states associated with strongly-coupled spins in 
spin clusters, as also reported in earlier studies >i2*2£ At 
still higher dilutions (not shown) , negative-energy states 



appear in the magnon spectrum, indicating instability 
of the collinear ferromagnetic state. The correlation be- 
tween spectral and spatial features of the magnon states 
is contained in the PR plots shown in Figure 3 for the two 
doping levels. In addition to the magnon softening with 
increasing dilution, the PR plots also show the increasing 
localization of magnon states with dilution, especially for 
the high-energy modes which are more likely to be local- 
ized in the strongly-coupled spin clusters. 

We note here that the magnon energy scale for ex- 
tended modes, which essentially determines the spin dy- 
namics in the ordered state due to thermal excitations 
of magnons, is very low compared to the hopping energy 
scale, as low as around ~ 0.02 for the x as 20% case in 
Fig. 2. On the other hand, corresponding to reduced 
fermion polarization in impurity-poor regions, the low- 
energy MF scale is J(aj)S/2 ~ 0.25. The nearly order- 
of-magnitude separation between these two energy scales 
implies, as discussed above, that the T = calculations 
provide a good description of the low-temperature MF 
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FIG. 5: Variation of T c with hole concentration p for a N — 
10 3 system with J = 4, x — 20%, and carrier bandwidth 
W = 12t — 10 eV, averaged over 50 configurations. 



state. 

Figure 4 shows effective control of impurity disorder 
at fixed dilution and doping. Enhanced potential disor- 
der with increasing J results in magnon softening (upper 
panel) , whereas negative Cd effectively reduces disorder in 
the doped spin-j band and results in magnon stiffening 
(lower panel), indicating stabilization of the ferromag- 
netic state and higher T c . 

Figure 5 shows the strong optimization of the tran- 
sition temperature T c with carrier concentration due to 
the characteristic competition between increasing overall 
magnitude of the carrier-induced spin couplings «/ 2 X?j 
and the increasing rapidity of its oscillation. Here we 
have estimated T c from the configuration average 




and taken a hole bandwidth (W = 12t = 10 eV) of the or- 
der of that of GaAs. Also, we have taken x = 20% for the 
impurity concentration (per unit cell for our sc lattice), 
corresponding to 5% Mn in the fee system Gai_ x Mn x As, 
which has 4 Ga sites per unit cell, and therefore impurity 
concentration Ax per unit cell. The peak at p/x ~ 1/5 
shows the importance of compensation in DMS systems, 
and the calculated T c is in good agreement with experi- 
mental results for Gai_ x Mn x As. 

The interplay of competing interactions and disorder 
is highlighted in Figure 6, which shows a comparison of 
magnon DOS and T c for the ordered and disordered cases 
at the same dilution x = 1/8. Here the ordered case cor- 
responds to a superlattice arrangement of impurities on 
alternate host lattice sites, studied earlier within a k- 
space sublattice-basis representation^ The low-energy 
part of the magnon spectrum, which corresponds to ex- 
tended modes, is seen to be significantly softened in the 
disordered case, resulting in substantially reduced T c , 
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FIG. 6: Comparison of magnon DOS and T c for the ordered 
and disordered cases with exactly same dilution x = 1/8, for 
a JV = 10 3 system with J = 4. Here p w 2.7% in the upper 
panel and W = 10 eV in the lower panel. 



and simply reflects the reduced spin stiffness due to com- 
peting interactions. The high-energy magnon modes in 
the disordered case are due to magnon modes localized 
over strongly-coupled impurity spin clusters, as already 
discussed above. The T c result for the ordered case is in 
close agreement with earlier fc-space analysis in terms of 
spin stiffnessr^ providing additional check on the validity 
of finite-size T c calculations. 



IV. FINITE-TEMPERATURE SPIN DYNAMICS 

In section III we saw that dilution-induced disorder re- 
sults in a strong enhancement in the density of low-energy 
magnons, with an appreciable fraction of localized modes 
corresponding to weakly-coupled spins, as well as forma- 
tion of impurity-spin clusters supporting localized high- 
energy magnon modes. In order to quantitatively inves- 
tigate the effect of these magnon features on the finite- 
temperature spin dynamics, in this section we evaluate 
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FIG. 7: Temperature dependence of averaged magnetization 
for different dilutions (x ~ 20%, 35%, 70%) and system sizes 
(N = 12 3 ,10 3 ,8 3 ) while keeping the number of magnetic 
states fixed. Here J = 4, W = lOeV, and p/x = 1/5. Also 
shown is the conventional spin- wave-theory result (25). 



the thermal reduction in magnetization due to magnon 
excitations. 

We first discuss a locally self-consistent magnon renor- 
malization scheme for the determination of local magne- 
tizations, which is equivalent to the site-dependent Tyab- 
likov decoupling procedure (local RPA) for a Heisenberg 
ferromagnet. We obtain the local magnetization from the 
Callcn formula 

, w _ (s g/Ki + $/) 25+1 + (s + 1 + $ 7 )$f +1 

W (l + ^+i.^s+i ( 22 ) 

for a quantum spin- S ferromagnet^ where we have in- 
troduced site-dependent boson occupation numbers 



1 



E 



(23) 



which explicitly involve the boson density {(j)\) 2 ■ A 
locally self-consistent magnon rcnormalization scheme 
(in terms of the fixed FKLM spin couplings Ju ~ 
(J 2 /4)[x°]/j) is then obtained if the boson Hamiltonian 
matrix elements (10-13) are self-consistently renormal- 
ized 



Tin 



E 



lx }u 2(55) 



(24) 



in terms of the local magnetization (Sf) instead of 
the MF values. Together with 7i\4>i) = u>i\<f>i), the 
coupled equations (22-24) then self-consistently yield 
the local magnetization (Sf) for all sites. In the 
translationally-symmetric case this yields the usual 



FIG. 8: Distribution of local magnetization, showing qual- 
itatively different spin dynamics for weakly- and strongly- 
coupled spins in a N — 10 3 system with J = 4, x — 20%, 
p « 4%, and W = lOeV. 



momentum-independent magnon-energy renormalization 
c ; = ({S')/S)uf. 

In the low-temperature regime, thermal magnon renor- 
malization is negligible ((Sf) — > S), and the bare (T = 0) 
magnon energies lj® provide a good description of the 
spin dynamics. We also have $/ << 1, and the magne- 
tization equation (22) reduces to the conventional spin- 
wave-theory result for the site-averaged magnetization 



(S z ) = S 



du> 



N(uo) 



(25) 



in terms of the bare magnon density of states N(uj). On 
the other hand, in an intermediate-temperature regime 
where k^T 3> for low-energy localized modes, then 
~ ksT/aji >> 1 and a nearly paramagnetic contri- 
bution (Sf) ~ uii/k-gT is obtained to the local magne- 
tization, highlighting the qualitatively different spin dy- 
namics of the isolated, weakly-coupled spins. 

Figure 7 shows the temperature-dependence of mag- 
netization evaluated from (22) and (23) with the bare 
(T = 0) magnon energies and wavefunctions and aver- 
aged over all impurity sites and several (20) configura- 
tions, for different dilutions (x w 20%, 35%, 70%) and 
fixed value of p/x = 1/5. Different host system sizes 
N = L 3 with L = 12, 10, 8 are taken so that the num- 
ber of magnetic states per configuration remains same 
(N m = 350) in Eq. (23). The strong enhancement in 
the thermal decay of magnetization with dilution reflects 
the effect of the large enhancement in density of low- 
energy magnetic excitations on the finite-temperature 
spin dynamics. Also shown (for x ps 20%) is the conven- 
tional spin- wave-theory result (25), which asymptotically 
approaches the previous result in the low-temperature 
regime, as expected. 

While the distinctly concave behaviour at higher tem- 
perature in the diluted case (Fig. 7, x = 20%) is in 
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FIG. 9: Configuration-averaged magnetization within an ap- 
proximate renormalization scheme involving the average bo- 
son occupation number, for a N — 12 3 system with J = 4, 
W = lOeV, x w 20%, and p « 4%, along with a fit with the 
Bloch form 5(1 - BT 3/2 ), with B = 0.6 x 10~ 3 KT 2/3 . 

itself not conclusive evidence for nearly paramagnetic 
behaviour of weakly-coupled spins, as similar behaviour 
would be obtained even for a translationally-symmctric 
system in a high-temperature (kftT 3> W; ) regime, the 
distribution of local magnetization (Sf) for a single con- 
figuration, clearly shows (Fig. 8) the qualitatively differ- 
ent spin dynamics for weakly- and strongly-coupled spins 
for the same magnon spectrum, thus highlighting the role 
of the spatial character of magnon states. The wide range 
of behaviour for different spins, and therefore the differ- 
ent local environments reflects the degree of complexity 
in a fully self-consistent magnon renormalization theory. 

We next include magnon renormalization within an 
approximate scheme^ suitable for the low-temperature 
regime where (Sf) ~ 5 for all sites. Here the average 
magnetization (S z ) is self-consistcntly calculated from an 
equation similar to (22) in terms of the configuration- and 
site-averaged boson occupation number 



* = (Si) = 




where the bare (T = 0) magnon energies uj? are uni- 
formly renormalized by the factor (S z )/S, as in the stan- 
dard Tyablikov theory. Figure 9 shows the temperature 
dependence of average magnetization obtained for sev- 
eral (N c = 20) configurations of a TV = 12 3 system with 
J = 4, x w 20%, and p w 4%. As expected, in the low- 
temperature regime where (S z ) — > 5 and << 1, the 
result approaches the previous (unrenormalized) result. 
Interestingly, the low-temperature behaviour of magneti- 
zation fits well with the Bloch form 5(1 - BT 3/2 ), with 
B = 0.6 x 10 -3 K~ 2 / 3 , which is of same order of mag- 
nitude (1.4xl0~ 3 K~ 2 / 3 ) as obtained in squid magne- 
tization measurements of an annealed 50nm sample of 
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FIG. 10: Substantial enhancement in T c due to spin clustering 
when the site-dependence of local magnon occupation number 
and magnetization is included, for a N — 10 3 system (single 
configuration), with J = 4, W = lOeV, x = 20%, and p w 4%. 

Gai_ x Mn x As with 6% Mn concentration! 2 ^. 

Dealing only with the average magnon occupation 
number, the above approximation does not incorporate 
the spatial segregation of magnon modes, particularly 
of the localization of high-energy modes over strongly- 
coupled impurity-spin clusters. Whereas, strong local 
correlations and high disordering temperature in these 
spin clusters should generally enhance long-range order- 
ing and T c , and therefore a quantitative analysis of this 
spin clustering is of interest. Localization of high-energy 
magnon modes over impurity clusters implies that clus- 
ter spins have relatively smaller participation in the low- 
energy modes, resulting in smaller magnon occupation 
number and higher local magnetization (Sf) on clus- 
ter sites at low temperature. Fig. 10 shows the substan- 
tial enhancement in T c resulting from including the site- 
dependence of S/ and (Sf), evaluated self-consistently 
from (22) and (23), but with a uniform magnon-energy 
renormalization loi = ((5 z )/5)w ; in (23) in terms of 
the average magnetization (S z ) = (l/N m )Y^i(3f)- Also 
shown is the result of the fully self-consistent magnon 
renormalization scheme (22-24), where the renormalized 
magnon energies loi and wave functions <f>i are obtained 
by diagonalizing the renormalized magnon Hamiltonian 
(24) which is updated at every step. 

V. FKLM WITH WEAKLY-COUPLED SPINS 

Even at x ks 20% impurity concentration, the domi- 
nant low-temperature magnetization behaviour (Fig. 9) 
in the diluted FKLM is determined by the majority low- 
energy extended modes, although a finite fraction of 
weakly-coupled spins is also clearly present (Fig. 8). In 
order to highlight the contribution of the low-energy lo- 
calized modes, in this section we consider a FKLM with a 
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FIG. 11: Magnon density of states and PR for the FKLM 
with a small fraction (x' = N' /N ~ 20%) of weakly-coupled 
spins (t' = 0.4) in a N = 8 3 system with J = 4 and p « 11%. 
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FIG. 12: Comparison of temperature dependence of local 
magnetizations (Sf) of a weakly-coupled spin (t' = 0.4) and 
a bulk spin in a TV = 8 3 system with J = 4 and p w 11%. 
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small fraction of weakly-coupled spins (resembling nearly 
paramagnetic impurities), which has some resemblance 
to the impurity-band model^S and evaluate the lattice- 
averaged magnetization within the locally self-consistent 
magnon renormalization scheme discussed in section IV. 
We therefore consider a FKLM 

I,S,<T " I 

where the nearest-neighbour hopping terms ti = t' < t 
for a small fraction (x' = N'/N) of randomly chosen 
sites, representing relatively isolated and weakly-coupled 
spins, within an impurity-band representation^ 

The magnon DOS and IPR (Fig. 11) clearly show low- 
energy localized modes corresponding to weakly-coupled 
spins, qualitatively similar to the diluted KLM case with 
randomly placed impurity spins, but with much sharper 
differentiation between extended (bulk) and localized 
(impurity) magnon states. The disorder here is rela- 
tively weaker and apart from the small fraction of low- 
energy magnon modes, remaining modes are essentially 
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FIG. 13: Lattice-averaged magnetization obtained within the 
locally self-consistent magnon renormalization scheme for the 
FKLM with N' = 150 weakly-coupled spins (t' = 0.3) in a 
N = 8 3 spin system with J = 4 and p ~ 11%. 



unaffected. For a single weakly-coupled spin in the lat- 
tice, Fig. 12 shows the qualitatively different (unrenor- 
malized) spin dynamics for the weakly-coupled and bulk 
spins obtained from Eqs. (22) and (23) with the bare 
magnon energies Results of the locally self-consistent 
magnon renormalization scheme (22-24) are shown in 
Fig. 13, which highlights both the dominant nearly para- 
magnetic spin dynamics of weakly-coupled spins in the 
low-temperature regime, and the dominant spin dynam- 
ics of bulk spins near T c . 
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FIG. 14: Magnon density of states for a N = 8 s system with 
different potential disorder strength W . Here J = 4 and hole 
doping p w 11%. 



VI. FKLM WITH POTENTIAL DISORDER 

Competing interactions between the randomly sepa- 
rated impurities due to positional disorder is an impor- 
tant feature of carrier-induced ferromagnetism in the di- 
luted FKLM. In order to highlight the interplay of po- 
sitional disorder and competing interactions, and distin- 
guish the role of purely potential disorder, we consider 
the concentrated FKLM with potential disorder 

H = t^2 a\ a a I+s ,a + e i a \,<j a h° - \ S i- (T i ( 28 ) 

7,5,(7 I,a I 

where the random on-site energies are distributed in the 
range — W/2 < ei < W/2. with a rectangular probabil- 
ity distribution of width W. Magnon energies are again 
obtained following the general method discussed in sec- 
tion II, by diagonalizing the magnon Hamiltonian (10) 
involving the effective spin couplings J 2 xu m terms of 
the particle-hole propagator xu- As there is no dilution, 
the number of magnetic states N m is equal to the num- 
ber of sites N. Without disorder, magnons in the con- 
centrated FKLM have been studied earlier in the context 
of heavy fermion materials^ ferromagnetic metals Gd, 
Tb, Dy, doped Eu& and manganitesi 29 i 30 i 31 i 32 

We find that purely potential disorder has a relatively 
milder effect on the magnon spectrum (Fig. 14), with 
a gradual shift in the magnon density of states peak to 
lower energy, indicating gradual decrease of T c with dis- 
order. While the magnon spectrum is substantially soft- 
ened, there are no negative-energy magnon modes even 
for disorder strength comparable to the bandwidth, and 
the collinear ferromagnetic state remains robust. The 
narrowing of the magnon energy distribution in Fig. 14 
is a typical feature of the disordered ferromagnet in which 
the low-energy extended modes are relatively less affected 
by the small averaged-out disorder on long length scales. 



By comparison, the interplay of competing interactions 
and disorder in the diluted case results in much stronger 
magnon softening (Fig. 2), and even with a comparable 
potential disorder ~ JS/2 = 5 for J = 4, increasing di- 
lution results in strong magnon softening and emergence 
of finite fraction of negative-energy modes indicating in- 
stability of the collinear ferromagnetic state. 

For the concentrated FKLM without disorder, the 
spin-stiffness expression consists of two parts — a posi- 
tive delocalization-encrgy term of order t favouring ferro- 
magnetic ordering and a negative exchange-energy term 
of order t 2 / JS favouring antiferromagnetic ordering. The 
magnon softening seen in Fig. 14 due to disorder can 
be readily understood in terms of a disorder-induced 
enhancement of the average exchange- energy term ~ 
(t 2 /(JS + e/ — e,/)), as also obtained for the disordered 
Hubbard model^i Essentially, enhancement in the ex- 
change term when the potential barrier is reduced (neg- 
ative e/ — e,/) more than compensates for its reduction 
when the barrier is enhanced (positive ej — ej), resulting 
in enhanced average exchange energy. 



VII. CONCLUSIONS 

The exact treatment of disorder in our finite-size cal- 
culations allowed for a quantitative investigation of the 
interplay of disorder and competing interactions in the 
carrier-induced ferromagnetic state of the diluted fer- 
romagnetic Kondo lattice model. Quantitative mea- 
sures of the competition between impurity spin cou- 
plings, stability of the ferromagnetic state, and the mag- 
netic transition temperature could be obtained from our 
real-space analysis of magnon excitations, with the esti- 
mated T c in good agreement with experimental results for 
Gai_ x Mn x As for corresponding impurity concentration, 
hole bandwidth, and compensation. Finite-temperature 
spin dynamics was studied within a locally self-consistent 
magnon rcnormalization scheme, equivalent to a site- 
dependent Tyablikov decoupling (local RPA), to investi- 
gate the role of magnon localization over isolated, weakly- 
coupled spins and over strongly-coupled spin clusters. 

The strong magnon softening observed with dilu- 
tion reflects increasingly effective competition between 
carrier-induced spin couplings, especially at lower hole 
doping where the couplings are longer-ranged. Indeed, 
the instability of the collinear ferromagnetic state at 
high dilution, as signalled by appearance of negative- 
energy magnon modes, indicates strong competition be- 
tween ferromagnetic and antiferromagnetic spin cou- 
plings. Highlighting the importance of compensation in 
DMS systems, the strong optimization of T c with hole 
doping at p/x << 1 stems from the competition between 
increasing overall magnitude and increasing rapidity of 
oscillation of the carrier-induced spin couplings. Com- 
parison of the ordered and disordered cases for exactly 
same dilution, with significant magnon softening and low- 
ering of T c in the disordered case, explicitly demonstrates 
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the interplay of disorder and competing interactions in 
DMS systems. The relatively milder magnon softening 
and robustness of ferromagnetic state obtained for the 
concentrated FKLM with purely potential disorder serves 
to highlight the greater degree of frustration due to po- 
sitional disorder in diluted systems. 

The large enhancement in density of low-energy mag- 
netic excitations due to competing interactions and dis- 
order is responsible for the strong thermal decay of mag- 
netization at high dilution. While the distribution of 
local magnetization clearly exhibits spin-dynamics be- 
haviour of both weakly- and strongly-coupled spins, the 
low-temperature behaviour of lattice-averaged magneti- 
zation was found to be dominated by low-energy ex- 
tended magnon states, and fitted well with the Bloch 
form Mo(l — BT 3 ' 2 ), with B of the same order of magni- 
tude as obtained in recent squid magnetization measure- 
ments on Gai_ x Mn x As samples^ 

The enhancement in ferromagnetism due to strong lo- 
cal correlations in strongly-coupled spin clusters is cap- 
tured in the locally self-consistent magnon renormaliza- 
tion scheme, as reflected in the substantial enhance- 
ment in T c obtained on including the site dependence 
of magnon occupation numbers. Furthermore, applied 
to the FKLM with a small fraction of weakly-coupled 
spins, the scheme highlighted both the dominant nearly- 
paramagnetic spin dynamics of weakly-coupled spins in 
the low-temperature regime, and the dominant spin dy- 
namics of bulk spins near T c . Thus a wide range of 
spin-dynamics behaviour involving different local envi- 
ronments can be studied practically within the self- 
consistent magnon renormalization scheme. 

Finally, we mention the outstanding and interesting 
issue of fermion-sector renormalization, in particular of 
the carrier-induced spin couplings due to quantum and 
thermal corrections to the bare particle-hole propaga- 
tor [x°(w)]. In the VCA-type approaches, the effec- 



tive impurity field seen by fermions is taken to vanish 
as T — > T c , yielding RKKY-type spin couplings. How- 
ever, the large separation between moment-melting and 
moment-disordering temperatures implies presence of ap- 
preciable local moments even near T c , so that fermions 
should continue to see impurity fields due to the slowly- 
fluctuating, locally-ordered impurity moments. Indeed, 
the splitting of fermion bands near T c even for rather 
moderate couplings (J/W ~ 0.2), obtained within dy- 
namical self-energy studies of the concentrated FKLM^ 
is precisely due to presence of the slowly-fluctuating im- 
purity fields. It will be of interest to examine both quan- 
tum and thermal corrections to spin couplings due to self- 
energy and vertex corrections in within a spin- 
rotationally invariant scheme which preserves the Gold- 
stone mode. Studied recently for the ferromagnetic state 
of the Hubbard modelf^l the net quantum correction can 
be essentially understood in terms of an exchange-energy 
correction due to fermion spectral- weight transfer. Pre- 
liminary calculation of spectral-weight transfer for the 
FKLM indicates a suppresion by the factor 1/5, sug- 
gesting relatively smaller quantum corrections for large 
spin quantum number 5. 
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